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ABSTRACT 

The recent years have seen combined measurements of X-ray and (weak) lensing contours for col- 
Uding galaxy clusters such as, for instance, the famous "Bullet" cluster. These observations have 
revealed offsets in the peaks of the baryonic and (dominant) gravitational matter component of order 
~ 100 — 200 kpc. Such discrepancies are difficult to explain using modified theories for gravity other 
than dark matter. Or are they not? Here we use the concept of "phantom dark matter" that is 
based upon a Newtonian interpretation of the MONDian gravitational potential. We show that this 
idea is in fact capable of producing substantial offsets in idealistic density configurations, involving 
a uniform external field. However, when analysed in a MONDian cosmological framework we deduce 
that the size (and probablity) of the effect is too small to explain the observed offsets found in the 
most recent observations, at least in the simplest incarnation of phantom dark matter as discussed 
here. The lensing centers in merging galaxy clusters are likely very close to the centers of true mass 
even in a MONDian cosmology. This gives the support to the idea that neutrino-like non-coUisional 
matter might be responsible for the observed offsets of lensing and X-ray peaks. 
Subject headings: galaxies: evolution - galaxies: halos - cosmology: theory - cosmology: dark matter 
- methods: iV-body simulations 



1. INTRODUCTION 

T he standard A Cold Dark Matter (ACDM) model 
(cf.. iKomatsu et aLll2008[ ) explains the formation of cos- 
mological structure in the non-linear regime in a hier- 
archical way, i.e. big structures are not formed mono- 
lithically but by successive merging of small struc- 
tures (e.g.. iDavis et al.lll985f) . Recent cosmological sim- 
ulations also support this idea o f hierarchical struc- 
ture formation in MOND gravity (iLlinares et al.l 2008 ) 



One unavoidbale consequence of any high speed colli- 
sion of mass concentrations seems to be the decoupling 
or offsetting of the baryonic component from the dark 
component. Aside from the aforementioned Bullet clus- 
ter - whose offset has been measured to be approxi- 
mately ^ 100 kpc (e.g. | Clo we et ar''2006f) - more ex- 



amp les are given in 



Jee et al.l ( 2005a ) 



IdlOOSb), 



and lBrada6 et al.l ( 2008h . The latter authors actually 



(but see also the ana lytical models of I Sander"^ (j2008D 
and IZhao erall (|2008l )). The hierarchical merging sce- 
nario naturally promotes the picture that we should ob- 
serve collisions of (clusters of) galaxies. Observation- 
ally there is evidence that some of these impacts actu- 
ally occur with speeds that are not readily reproduced 
by simulations of ACDM structure formation, in the 
sense that the relative speed of the merging dark ha- 
los is rarely very higher than the internal dispersion 
of each halo dH avashi fc White 2006; 'Kncbe et alll2008t 
iLhnares et alB oOS: Angus & McCaugh 2008). There is, 
for instance, the famous "Bullet cluster" , an extremely 
high velocity merger between two galaxy clusters, with 
an inferred shock velocity of ~ 4700km/sec. While rela- 
tive encounters with comparable velocities are rather rare 
in pure dark matter sim ulations (e.g. Havashi & Whit^ 
^2Xm . lKnebe"eraII l2008h . they may nevertheless be ac- 
comodated when considering explicit hydrodynamic al 
modelling of the phenomenon (|Springel fc Farrarll2007t ) . 
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present data for a particulr cluster (i.e. MACS J0025.4- 
1222) with an even greater difference of ~ 200 kpc in 
the peak of the baryonic and the dark matter. All this 
work is cu l minat ed in one particular recent Letter by 
IShan et ahl ()2009f ) where a sample of 38 galaxy clusters 
has been studied utilizing both X-ray and strong lensing 
observations. They show that such offsets are a com- 
mon phenomenon in galaxy clusters: they found at least 
13 objects with a separation greater than 50 kpc with 
3 clusters exhibiting a separation between baryonic and 
(hypothetic) dark component in excess of 200 kpc. 

All these papers are analysing combined X-ray and 
lensing observations to decipher (and actually measure) 
the offset between baryonic and dark matter. But how 
certain are we that the lensing signal is caused by "real" 
dark matter particles? What if the gravitational po- 
tential is not generated by Newtonian physics yet inter- 
preted in that way? 

One theory capable of producing potentials akin to 
dark matt er particles is modified Newtonian dynamics 
(MOND, lMilgrom|[l98l : (pedagogical) reviews of the 
concepts and successes can be found in , for i nstance, 
ISanders fc McGaughl (poM ) and lMilgromI (l2008h . While 
MOND was originally proposed as an alternative to New- 
tonian gravity designed to solely explain galactic dynam- 
ics without the need for dark matter the theory has 
gained substantial momentum during the past decade: 
although current cosmological observations point to the 
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existence of vast amoun ts of non-baryonic d ark mat- 
ter in the Universe (e.g. iKomatsu et al.l[2008f ). not all 
of the features of CDM models appear to match ob- 
servational data (e.g., the "missing satellite problem" 
(jKlvpin et al.lll999l : iMoore et all ' 19991 an d the so-called 
"cusp -core crisis" (e.g. lde Blok et al . 20031 ISwaters et"ari 
120031) ). Just as for CDM the MONO theory success- 
fully matches observations on a wide range of scales, 
different types of galaxies i ncluding dwarfs and gi- 
ants, spirals and elhpticals (iFamaev &: Bin ncv 2005'; 
Gentile et all 120071: iMilgrom fc Sanderd l2007c Milgrom 



200^ Sanders fc Noordermeeil 120071 : iMcGaughl I200E 



Angusli2008.) . However, one of MOND's major set-backs 



for a long time was the lack of a covariant formulation 
of the theory. This has been remedied by iBekensteinI 
((200l who was the first to cast MOND into a more 
universal form compliant with general relativity. This 
in turn spawned further investigations into the same 
direction leaving us nowadays with variou s relativis- 
tic fo rmulations of t he MOND theory (e.g. [B ekcnsteiij 
20041 'Sanders 20051; iBr uneton fc Esposito-Farese 2007i 



Zhao 2007; Zlosnik et al. 2007; Zhao '20081: ISkordidlloOS : 
Blanchet fc Le Tieo .20091) : a recent review of both 



MOND and its relati vistic offspr i ng (in particular the 
T eVeS for mulation of Bckcnstcin (2004!)) can be found 
in IMilgrom (2008) and Skordis (2009). We though need 
to acknowledge that despite the original idea of aban- 
doning the need for dark matter, even MOND cannot 
do without it: despite the great success we need to 
accept that even MOND cannot do well without dark 
matter completely. A recent study utilizing a combi- 
nation of strong and weak lensing by galaxy clusters 
indicates the necessity fo r neutrinos of mass 5 — 7eV 
(jNatarajan fc Zhaol l2008f l. And to be consistent with 
dark matter estimates of galaxy clusters and observa- 
tions of the CMB anisotropics Angus (2009); Angus et al. 
(|2009) claims for lleV^ neutrinos. One theory capabale of 
accomodating both t hese r e quirem ents is that of a mass- 
varying neutrino by IZhad (j2008l ). In summary, we are 
eventually left with a situation where the development 
of several frameworks for a relativistic formulation of 
MOND enabled the study of the cosm i c mic rowave back- 
ground (jSkordis et gd] 120061: ETet al.l 120081), cosmolog i- 
cal structure formation (jHalle et al.ll2008l : lSkordiil2008h. 
strong gra vitational lensing of galaxies ( Zhac l2006t 
IChen fc Zh ao 2006; Shan et al. 200 8), and weak l ensing 
of clusters (|Angus et all 120071 : IFamaev et al.ll2008f) . The 
MOND theory has matured and became a credible com- 
petitor to the commonly accepted CDM model. 

In that regards, however, it appears important to look 
for even more tests that are capable of discriminating 
between MOND and Newtonian gravity, especially in 
the context of cosmology. We therefore raise the ques- 
tion whether the kinds of offsets alluded to above can 
be explained by simply interpreting the MONDian po- 
tential in a Newtonian way? This concept of "phantom 
dark matter" was originally introduced in the beginnings 
of MOND alread y by Milg rom (1986 ) and has rec ently 
been discussed bylMilgroml (|2008). Wu et all (|2008[) . and 
iBienavme et al.l ()200 9^. It is based upon the idea to use 
the MONDian potential in a dark matter context, i.e. 
given the MONDian potential one can use the Newto- 
nian Poisson's equation to derive the corresponding den- 
sity of matter that would be needed in the Newtonian 



context. Then, subtracting the visible (baryonic) mat- 
ter one obtains the "virtual" dark matter or, in other 
words, "phantom dark matter" distribution predicted by 
MOND. And in a Newtonian interpretation this phantom 
dark matter would be responsible for the gravitational 
lensing signal alluded to above. 

We though n e ed to acknowledge that 
iBrownstein fc MoffatI (|2007f ) already pointed out the 
possibility that the observed offset in the (alleged) dark 
matter and baryonic density peaks of the Bullet cluster 
system can be explained by extending the equations for 
gravitational lensing to modified gravity, without the 
need for a dominant dark matter component. Further, 
as MOND is a non-linear theory it is not clear whether 
the (baryonic) matter will be distributed ab initio in 
the same way as phantom dark matter. We therefore 
set out to answer the question whether or not these 
two density fields share peaks at the same locations 
and are distributed in comparable ways, respectively. 
Can an offset between the dark and baryonic matter 
be explained by the non-linearity of the MONDian 
Poisson's equation and the existence of this putative 
phantom matter, respectively? 

We need to close with a cautionary note: this work 
does not deal with (collisions of) galaxies or galaxy clus- 
ters; we are solely focusing on the properties of the mat- 
ter density fields and their respective peaks. The primary 
question we set out to answer is whether or not MOND 
will produce offsets between the actual (baryonic) matter 
component and the phantom matter field, even though 
this work is motivated by observations of such offsets in 
galaxy clusters. 

2. THE NON-COSMOLOGICAL FRAMEWORK 

Before investigating phantom dark matter in a cos- 
mological environment we start off with phrasing the 
question about shifts in the respective density peaks for 
MONDian systems in a non-cosmological context. This 
will provide us with a gauge whether or not we should 
actually expect to find the reputed offsets. 

2.1. Phantom Dark Matter 

The MONDian Poisson's equation embedded within an 
external field reads as follows 



V 



ao 



g = gcxt - V<I>int, (1) 



where p is the baryonic matter, goxt an external field, 
and <I>int the (internal) potential of the system and p.{x) 
is /i — > 1 for a; ^ 1 (Newtonian limit) and p —>■ x for 
a; <C 1 (deep MOND hmit).^. We now take the liberty to 
interpret this internal potential $i,it within the context 
of Newtonian gravity 



V^^int =47rG(/9-|-pph), 



(2) 



where the (p + pph) is the total dynamical mass of the 
system, and pph is the so called "phantom dark mat- 
ter" , which can be held responsible for the "extra grav- 
ity beyond the baryonic matter"^ in the linear Poisson's 

^ Please n ote that we use d = x(l + x^)~^^^ as originally 

suggested by lMilgroml 1)19831 ) throughout our tests. 
^ Or in other words "dark matter" . 
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TABLE 1 

Parameters of our non-cosmological test model. The mass 

IS IN UNIT OF (IO^^Mq), positions AND Tq ARE IN UNITS OF kpc, 
AND THE EXTERNAL GRAVITY ACCELERATION IS IN UNIT OF aQ. 



Model 


Mass 


Centre 


9cxt 


TO 


G+EF 


10.0 


(0,0,0) 


1.0 


2.0 



equation Eq. But as opposed to the dark matter 

theory, MOND immediately predicts the distribution of 
the dynamical mass as soon as the baryons p and the 
gravity of the environment goxt are specified. However, 
due to the external field goxt the boundary condition of 
the (internal) system changes, not necessarily preserving 
spherical symmetry: the distribution of the dynamical 
mass is somewhat d ifferent from that of CDM halos (see 
IWu et al.ll2007l .'2008). Further, there exist negative solu- 
tions of the phantom dark matter, and the peaks of dy- 
namical mass can in fact be offset to the baryonic peaks! 
These effects are most significant at places where the ex- 
ternal and internal fields are comparable and should be 
quantified more carefully in the following subsection. 

2.2. The Simulations 

For the non-cosmological settings studied in this Sec- 
tion we use the MON Dian Poisson solver developed 
by th e Bologna group (jCiotti et al.l 120061 : iNipoti et al.l 
|2007| ) to solve Eq. ([T]) and hence derive the internal 
potential $int of the systems under investigation. The 
Poisson solver is a spherical grid code, and our choice 
for the grid parameters is rir x ng x — 256 x 
64 X 128 with a radial grid spacing given by r.; — 
ro tan [(i -I- 0.5)0.57r/(nr -I- 1)] kpc. We further utihze 
Eq. ^ to derive the dynamical mass and the phantom 
matter, respectively. 

All our (baryonic) galaxies have a Plummer density 
profile 



p{r) 



f 3M \ 
I 4^ j 



r 



-5/2 



(3) 



with a core radius b of 1.0 kpc. 



We investigated several scenaria (cf. discussion below 
in Section [231) , however, decided to only present the re- 
sults for one representative configuration where a single 
galaxy G is embedded within a strong constant exter- 
nal field EF. A summary of the actual parameters of the 
model can be found Table [TJ 

2.3. Density Peak Ojfsets 

Figure [T] shows the distribution of the dynamical mass- 
for model G-I-EF, i.e. a single galaxy embedded within 
an external field of strength |gcxt| = 0,0 along the x-axis. 
We notice that when the internal and external field are 
of the same order of magnitude, i.e. at x ~ 15 kpc there 
are some noticable effects: first, we obtain negative dy- 
namical mass on the corns perpendi cular to t h e dire ction 
of external field also mentioned by iMilgron] (|1986[ ) and 
iWu et al] ()2008f ): second, there exist an additional peak 
on the X-axis, right where the external and internal fields 
cancel each other! However, this additional peak is four 
orders of magnitudes smaller than the baryonic peak. 




-10 

X (kpc) 

Fig. 1. — Isodensity countours of the dynamical mass, i.e. 
baryons + phantom dark matter density, on the x — z plane for 
a galaxy embedded within an external field along the x-axis. 

We have also run idealistic simulations with (un-)equal 
mass galaxies with and without external fields (not 
shown though). In summary, we have seen that in most 
of the situations the interpretation of the MONDian po- 
tential in a Newtonian sense will lead to the prediction 
of additional peaks in the distribution of the dynamical 
mass when compared to the actual (underlying) baryonic 
matter distribution. However, the strength of these extra 
peaks vary and depend on the actual setup of the system 
ranging from as low as four orders of magnitude smaller 
to as large as 1% for the cases considered here. 

Encouraged by the observation that we actually re- 
cover offsets in our controlled experiments we may now 
rightfully ask the question whether these additional phan- 
tom peaks occur in realistic cosmological simulations^ i.e., 
whether a self-consistent cosmological simulation will 
provide a suitable variety of configurations so that we 
will in fact be able to observe (and quantify) the offset 
between the baryonic and phantom matter peaks. 

3. THE COSMOLOGICAL FRAMEWORK 

3.1. Phantom Dark Matter 

The equation for the (MONDian) gravitational poten- 
tial in a cosmological setting is somewhat different to 
Eq. ([T]) and reads 



|V$ 



07(0) 



M 



AttG 



ip-p), (4) 



We further took the liberty to encode the MONDian ac- 
celeration scale 7(a) as a (possible) function of the cos- 
mic expansion factor a. The most naive choice would 
be 7(0) = go = 1.2 x 10~*cm/sec^ whereas other theo- 
ries may lead to different dependencies; for instance, in 
IZhaol ()2008f ) 7(a) is given as 7(a) = a^^'^go. For more de- 
ta ils and a deriva t ion of this equation we refer the reader 
to iLlinares et al.l ()2008[ ) where it has been justified and 
implemented into t he cosmological A^-body code MLAPM 
(jKnebe et al.l l200n. 

Given the MONDian potential we may now ap- 
ply the same logic as in Section [2] and use it with the 
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Newtonian Poisson's equation whose right-hand-side will 
again no longer be the (baryonic) density p field alone but 
rather read as follows 



AttG 



[P + Pph) -{p + Pph) 



(5) 



This is the defining equation for the phantom matter 
density field pph used througout this Section. 

3.2. The Simulation 

The analysis presented in h ere is based upon a par- 
ticular simulation published in lLlinares et al.l (|2008D . i.e. 
the OCBMond2 model. This simulation has been run in 
a cosmological volume with a side length of 32/i~^Mpc 
and utilized 128'^ particles. It employe d the MONDi- 
fication of the TV-b ody code MLAPM (Kn ebe et al.l 120011 : 
iLHnares et al.l[200^ . We chose to simulate an open uni- 
verse with neither dark matter nor dark energy but char- 
acterized by fif, = 0.04. The simulation was started at 
redshift z = 50 and resorted to a Hubble parameter of 
h = 0.7. We further need to mention that there are 
two values tJg in a- MOND simulation, one characterising 
the amplitude of fluctuations of the initial condition and 
one measuring the strength of fluctuations at the present 
time. This comes about bec ause of the faster growth of 
struc tures in MOND (cf. [Sand ers 2001; Kneb(L& Gibsoii| 
I2OOI iLhnares et aP [2001 7 i.e. in order to arrive at a 
comparable evolutionary stage to a ACDM model at red- 
shift z = with (Tg ~ 0.9 we had to lower the magnitude 
of the fluctuations to a% = 0.4 during the process of gen- 
erating the initial condition. We acknowledge that this 
value is incompatible with CMB constraints, at least in 
the dark matter explanation for cosmic structure forma- 
tion. MOND, however, is a highly non-linear theory and 
the simulation presented and used here should be consid- 
ered as a first toy model for trying to understand struc- 
ture formation using modified gravity. For more details 
and elab orate study of t he sim ulation we refer the reader 
again to lLlinares et al.l (|2008f ). 

In order to perform the analysis presented here we fur- 
ther modified our potential solver to use the MONDian 
potential obtained by solving Eq. ((4]) for the final 
output at today's redshift z — Q and our knowledge about 
the (baryonic) matter density p together with Eq. (O to 
obtain the resulting phantom density pph. 

3.3. Locating Density Peaks 

Given both the (baryonic) matter density p and the 
phantom matter density pph we determine peaks in both 
fields. To this extend we smooth the fields on various 
scales and study spatial offsets in corresponding peaks 
in relation to the smoothing scale. This is accomplished 
by exploiting the adaptive mesh nature of the simulation 
code MLAPM used to generate the simulation in the first 
place (jKnebe et al.ll2001l : ILhnares et al.ll2008[ ): the (bary- 
onic) matter is represented by discrete particles whose 
mass is assigned to a regular grid covering the whole 
computational volume. This grid is then recursively re- 
fined in regions of high (baryonic) density according to 
a pre-selected refinement criterion of 4, 8, or 16 particles 
per cell. This leaves us with a hierarchy of (nested) re- 
finement patches where the boundary of each such patch 
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Fig. 2. — Sketch illustrating the definition of the (baryonic) 
matter density centres x\ on various refinement patches and the 
corresponding phantom matter peaks a;?*^. Note that the boundary 
of each patch is an isodensity contour p^"" in the (baryonic) matter 
distribution. Due to the nature of the mass assignment of the 
(baryonic) particles onto each refinement patch we are left with a 
density field smoothed on approximately the scale of the respective 
grid spacing ti. 



TABLE 2 
Smoothing scales in h~ 



[fe-ikpc] 



16382 
8192 
4096 
2048 
1024 
512 



1.95 
3.91 
7.81 
15.63 
31.25 
62.50 



defines a unique (baryonic) isodensity contour. The sit- 
uation is illustrated in Figure |2| where we show an exam- 
ple of two nested refinements embedded within a regular 
domain grid covering the whole computational volume 
of side length B. For each isolated patch we calculate 
x\ and by first finding the position of the cell con- 
taining the maximum in p (pph); we then use the density 
(phantom density) weighted average of the 27 neighbour- 
ing cells to define x\ {xf^). The relevant quantity to be 
studied below is the difference between the two peaks in 
the density field 



D 



X 



phi 



(6) 



that obviously is both a function of the smoothing scale 
€i and the isodensity contour pf°. 

3.4. Density Peak Offsets 

We are primarily interested in the question whether 
or not there are any (substantial) offsets in the peaks of 
the (baryonic) matter density field and the corresponding 
phantom field defined via Eq. ([5]). For the time being 
we therefore ignore any relation this offset has with the 
refinement level it is based upon (cf. Eq. (O). In Figure|3| 
we simply plot the cumulative distance distribution N{< 
D) normalized to the total number of refinement patches 
on all levels; we further chose to normalize the distance to 
the respective smoothing scale as we consider distances 
smaller than this scale to be below the resolution limit 
and hence not credible. The physical values for for 
the grid levels used in our calculations are summarized 
in Table H 
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Fig. 3. — The cumulative number distribution A'(< D) of the 
offset D between (baryonic) matter density and phantom density. 
The offset has been normalized to the respective smoothing scale 
of the refinement patch it is based upon. 




0.1 1.0 
distance [snnoottiing scale] 

Fig. 4. — The relation between the (baryonic) isodensity level of 
the refinement patch and the (normalized) offset between matter 
and phantom density on it. 

While most of the offsets between (baryonic) matter 
and phantom matter are in fact smaller than the resolu- 
tion limit there are nevertheless of order 1% instances for 
which we observe larger (and hence physical) differences! 

We like to caution the reader that the same matter 
peak enters multiple times (at most six times) into Fig- 
ure[3](and all subsequent plots below). This is due to the 
fact that we smooth the same peak using various smooth- 
ing scales ei listed in Table [2l However, as we are not 
interested in the change in offset for a given peak when 
altering the smoothing scale we can treat them indepen- 
dently. 

Even though we just found that there is a small yet 
measurable probability of finding an offset between bary- 
onic and phantom density it still remains unclear how 
this can be interpreted in terms of astrophysical objects. 
Observationally the edge of an object is primarily defined 
by a given threshold in (over-) density. This, however, is 
a natural by-product of our method for calculating D (cf. 



0.0100 





2048 


1024*\ 


1 

512 \ 











10 100 
distonce [kpc/li] 

Fig. 5. — Cumulative distribution of isolated patches with an 
offset D > (normalized to the total number of patches). Only 
patches on grids L < 2048 fulfill this criterion and the contribution 
of these grids to the distribution are marked by the respective L 
values given in the plot. 

Section 13. 4p . We therefore plot in Figure [4] the depen- 
dence of the (normalized) distance D on the overdensity 
of the corresponding refinement patch. Recall that the 
usual overdensity limit for virialized objects in an ACDM 
cosmology at redshift z = is w 340 and approximately 
coincides with our coarsest refinement level. 

As the only credible difference in the position of (bary- 
onic) matter and phantom matter peaks are apparent on 
the lower iso-density levels (i.e. the physically larger re- 
finement patches) one may raise the question whether we 
calculated the offset for corresponding peaks. A large re- 
finement patch will certainly host several peaks both in 
baryonic and phantom matter, so how can we be sure to 
take the difference between matching peaks? Maybe the 
maximum baryonic density is not at the same position as 
the maximum phantom density (cf. Section l3.4p ? This 
concern is readily eliminated as the maximum offset ob- 
served is no larger than two times the actual smoothing 
scale, i.e. the peaks lie in two neighbouring grid cells. 

So far, we always normalized the offset D to the re- 
spective smoothing scale ei. However, to gain a better 
and more quantitative feeling for the relevance of our re- 
sults it appears obligatory to also consider the distance in 
physical units /i^^kpc. To this extent we plot in Figure[5] 
the cumulative distribution of all offsets D > larger 
than the respective smoothing scale in physical h^^kpcJ 
Note that only the three coarsest grids (i.e. 512"^, 1024'^, 
and 2048"^) lead to offsets that are larger the smoothing 
scale; for all finer grids the distance in the (baryonic) 
matter and the phantom matter peak is below the credi- 
bility level given by the smoothing scale. We further ob- 
serve that the reasonable offsets lie in the range between 
15-80/i~^kpc. We though need to acknowledge that the 
absolute fraction of isolated patches fulfilling this credi- 
bility criterion is below 0.6%. Our simulation therefore 
has difficulties to accomodate offsets as large as the ones 
observed. 

^ Note that this figure represents a zoom of Figure \3\ into the 
region right of the vertical line and in physical units on the x-axis 
now. 
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4. SUMMARY AND CONCLUSIONS 

Driven by the observation of offsets in the baryonic and 
gravitational m atter distribution in coUisions of galaxy 
clusters (e.g. iShan et all l2009f) we explore such phe- 
nomena in the context of phantom dark matter. This 
is an interpretation of the MONDian potential (gener- 
ated purely by baryons) in a Newtonian context, i.e. 
the MOND potential is used with the standard (New- 
tonian) Poisson's equation and the resulting right-hand- 
side source term is understood as a combination of the 
baryonic matter and some phantom dark matter. 

An initial study of (interacting) galaxies in isolation as 
well as external fields indicated that we should expect to 
find additional peaks in the distribution of the dynam- 
ical mass as opposed to the baryonic mass distriution. 
However, the strength of these extra peaks varied and 
the contrast may in some cases be even too low to be 
observed. 

Utihzing a MON Dification of the TV-body code MLAPM 
(|Knebe et al.ll200l[) we set ourselves into the position of 
calculating both the baryonic and phantom density dis- 
tributions in a fully self-consistent MONDian cosmolog- 
ical simulation on adaptive refinement patches. We then 
quantified differences in the peaks of both these fields 
concluding that the (theoretically predicted) offsets are 
too small to be compliant with the observed offsets, at 
least in the presented incarnation of phantom matter and 
our MONDian cosmological simulation. 



One possible drawback of our applied method is the 
fact that the isodensity level that define isolated refine- 
ment patches are based upon pf, only. However, we com- 
pensated for this quibble by adjusting the refinement cri- 
terion and subsequentially modifying the size of isolated 
patches; we though could not detect any systematics. 

We conclude that our results give support to the idea 
that neutrino-like non-coUisional matter might be re- 
sponsible for the observed offsets of lensing and X-ray 
peaks. There are in fact indications by several au- 
thors that non-classical neutrinos are require d to explain 
phenomena, such as, cluster l ensing ^N atarai an fc Zhaol 
I2008D or CMB anisotropics ()Angus et all i200i ) within 
the context of (relativistic) MOND. One theory capa- 
bale of accomodating both these requirements is that of 
a mass- varying neutrino by Zhao (.20081 ) to be studied 
more detailed in future work. 



AK is supported by the MICINN through the Ra- 
mon y Cajal programme. CL and AK further acknowl- 
edge funding by the DFG under grant KN 755/2. This 
work was also carried out under the HPC-EUROPA-f-|- 
project (project number: 211437), with the support of 
the European Community - Research Infrastructure Ac- 
tion of the FP7 "Coordination and support action" Pro- 
gramme. XW thanks the support of the SUPA stu- 
dentship. 



REFERENCES 



Angus, G. W. 2008, MNRAS, 387, 1481 
— . 2009, MNRAS, 394, 527 

Angus, G. W., Famaey, B., & Diaferio, A. 2009, ArXiv e-prints 
Angus, G. W., & McGaugh, S. S. 2008, MNRAS, 383, 417 
Angus, G. W., Shan, H. Y., Zhao, H. S., & Famaey, B. 2007, ApJ, 
654, L13 

Bekenstein, J. D. 2004, Phys. Rev. D, 70, 083509 
Bienayme, O., Famaey, B., Wu, X., Zhao, H. S., & Aubcrt, D. 2009, 
ArXiv e-prints 

Blanche!, L., & Le Tiec, A. 2009, Phys. Rev. D, 80, 023524 
Bradac, M., Allen, S. W., Treu, T., Ebeling, H., Massey, R., Morris, 

R. G., von der Linden, A., & Applegate, D. 2008, ApJ, 687, 959 
Brownstein, J. R., & Moffat, J. W. 2007, MNRAS, 382, 29 
Bruneton, J.-P., & Esposito-Farese, G. 2007, Phys. Rev. D, 76, 

124012 

Chen, D.-M., & Zhao, H. 2006, ApJ, 650, L9 

Ciotti, L., Londrillo, P., & Nipoti, C. 2006, ApJ, 640, 741 

Clowe, D., Bradac, M., Gonzalez, A. H., Markevitch, M., Randall, 

S. W., Jones, C, & Zaritsky, D. 2006, ApJ, 648, L109 
Davis, M., Efstathiou, G., Frenk, C. S., & White, S. D. M. 1985, 

ApJ, 292, 371 

de Blok, W. J. G., Bosma, A., & McGaugh, S. 2003, MNRAS, 340, 
657 

Famaey, B., Angus, G. W., Gentile, G., Shan, H. Y., & Zhao, 
H. S. 2008, in Proceed, of the Sixth International Heidelberg 
Conference on Dark Matter in Astro- and Particle Physics 
(DARK 2007), held at the University of Sydney, Australia, 23-28 
September 2007, Eds. H. V. Klapdor-Kleingrothaus, G. Lewis, 
World Scientific, Singapore, ISBN-13 978-981-281-434-0, ISBN- 
10 981-281-434-5, pp. 393-401, 393-401 
Famaey, B., & Binney, J. 2005, MNRAS, 363, 603 
Gentile, G., Famaey, B., Combes, F., Kroupa, P., Zhao, H. S., & 

Tiret, O. 2007, A&A, 472, L25 
Halle, A., Zhao, H., & Li, B. 2008, ApJS, 177, 1 
Hayashi, E., & White, S. D. M. 2006, MNRAS, 370, L38 
Jee, M. J., White, R. L., Benftez, N., Ford, H. C, Blakeslee, J. P., 
Rosati, P., Demarco, R., & lUingworth, G. D. 2005a, ApJ, 618, 
46 

Jee, M. J., White, R. L., Ford, H. C, Blakeslee, J. P., lUingworth, 
G. D., Coe, D. A., & Tran, K.-V. H. 2005b, ApJ, 634, 813 



Klypin, A., Gottlober, S., Kravtsov, A. V., & Khokhlov, A. M. 

1999, ApJ, 516, 530 
Knebe, A., Arnold, B., Power, C, & Gibson, B. K. 2008, MNRAS, 

386, 1029 

Knebe, A., & Gibson, B. K. 2004, MNRAS, 347, 1055 
Knebe, A., Green, A., & Binney, J. 2001, MNRAS, 325, 845 
Komatsu, E., Dunkley, J., Nolta, M. R., Bennett, C. L., Gold, 
B., Hinshaw, G., Jarosik, N., Larson, D., Limon, M., Page, L., 
Spergel, D. N., Halpern, M., Hill, R. S., Kogut, A., Meyer, S. S., 
Tucker, G. S., Weiland, J. L., WoUack, E., & Wright, E. L. 2008, 
ArXiv e-prints, 803 
Li, B., Barrow, J. D., Mota, D. F., & Zhao, H. 2008, Phys. Rev. D, 
78, 064021 

LUnares, C, Knebe, A., & Zhao, H. 2008, MNRAS, 391, 1778 

McGaugh, S. S. 2008, ApJ, 683, 137 

Milgrom, M. 1983, ApJ, 270, 365 

— . 1986, ApJ, 306, 9 

— . 2007, ApJ, 667, L45 

— . 2008, ArXiv e-prints 

Milgrom, M., & Sanders, R. H. 2007, ApJ, 658, L17 

Moore, B., Ghigna, S., Governato, F., Lake, G., Quinn, T., Stadel, 

J., & Tozzi, P. 1999, ApJ, 524, L19 
Natarajan, P., & Zhao, H. 2008, MNRAS, 389, 250 
Nipoti, C, Londrillo, P., & Ciotti, L. 2007, ApJ, 660, 256 
Sanders, R. H. 2001, ApJ, 560, 1 
— . 2005, MNRAS, 363, 459 
— . 2008, MNRAS, 463 

Sanders, R. H., & McGaugh, S. S. 2002, ARA&A, 40, 263 
Sanders, R. H., & Noordermeer, E. 2007, MNRAS, 379, 702 
Shan, H. Y., Feix, M., Famaey, B., & Zhao, H. 2008, MNRAS, 387, 
1303 

Shan, H.-Y., Qin, B., Fort, B., Tao, C, & Wu, X.-P. 2009, ApJ, 
submitted 

Skordis, C. 2008, Phys. Rev. D, 77, 123502 

— . 2009, Classical and Quantum Gravity 26, 143001 

Skordis, C, Mota, D. F., Ferreira, P. G., & Boehm, C. 2006, 

Physical Review Letters, 96, 011301 
Springel, V., & Farrar, G. R. 2007, MNRAS, 380, 911 
Swaters, R. A., Madore, B. F., van den Bosch, F. C, & Balcells, 

M. 2003, ApJ, 583, 732 



7 



Wu, X., Famaey, B., Gentile, G., Perets, H., & Zhao, H. 2008, 

MNRAS, 386, 2199 
Wu, X., Zhao, H., Famaey, B., Gentile, G., Tiret, O., Combes, P., 

Angus, G. W., & Robin, A. C. 2007, ApJ, 665, LlOl 
Zhao, H. 2006, ArXiv Astrophysics e-prints 
— . 2007, ApJ, 671, LI 



— . 2008, ArXiv e-prints, 805 

Zhao, H., Xu, B.-X., & Dobbs, C. 2008, ApJ, 686, 1019 
Zlosnik, T. G., Ferreira, R G., & Starkman, G. D. 2007, 
Phys. Rev. D, 75, 044017 



